o 
o 

(N 



> 



o 



X 



On the abundance discrepancy problem in H II regions. 1 

Jorge Garcia- Roj as 
Institute) de Astrofisica de Canarias, E-38200, La Laguna, Tenerife, Spain 

j ogarcia@ll . iac . es 
and 
Cesar Esteban 

Instituto de Astrofisica de Canarias, E-38200, La Laguna, Tenerife, Spain 

cel@ll . iac . es 

ABSTRACT 



Si 

6 

The origin of the abundance discrepancy -i.e. the fact that abundances de- 



rived from recombination lines are larger than those from collisionaly excited 
lines- is one of the key problems in the physics of photoionized nebula. In this 
work, we analize and discuss data for a sample of Galactic and extragalactic H II 
i/-) ! regions where this abundance discrepancy has been determined. We find that 

the abundance discrepancy factor (ADF) is fairly constant and of the order of 2 
in all the available sample of H II regions. This is a rather different behaviour 



than that observed in planetary nebulae, where the ADF shows a much wider 
range of values. We do not find correlations between the ADF and the O/H, 



0++/H+ ratios, the ionization degree, T e (High), T e (Low)/ T e (High), FWHM, 
and the effective temperature of the main ionizing stars within the observational 
uncertainties. These results indicate that whatever mechanism is producing the 
abundance discrepancy in H II regions it does not substantially depend on those 
nebular parameters. On the contrary, the ADF seems to be slightly dependent 
on the excitation energy, a fact that is consistent with the predictions of the 
classical temperature fluctuations paradigm. Finally, we obtain that T e values 
obtained from O II recombination lines in H II regions are in agreement with 
those obtained from collisionally excited line ratios, a behaviour that is again 
different from that observed in planetary nebulae. These similar temperature 



1 Based on observations collected at the European Southern Observatory, Chile, proposals number ESO 
68.C-0149(A) and ESO 70.C-0008(A). 
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determinations are in contradiction with the predictions of the model based on 
the presence of chemically inhomogeneous clumps but are consistent with the 
temperature fluctuations paradigm. We conclude that all the indications suggest 
that the physical mechanism responsible of the abundance discrepancy in H II 
regions and planetary nebulae are different. 

Subject headings: line: identification. ISM:abundances — H II regions. 



1. Introduction 

H II regions are essential tools for the study of the chemical composition and star for- 
mation in the Universe, specially in the extragalactic domain. It is necessary to be confident 
that our traditional techniques for deriving chemical abundances in ionized nebulae (based 
on the analysis of intensity ratios of colisionally excited lines, CELs) provide the real values. 



Forty years ago, Peimbert ( |1967| ) characterized the temperature structure of a gaseous 



nebula to a second order by two parameters: the average temperature, T and the mean 
square temperature fluctuation, t 2 . At least, two independent methods of determining the 
electron temperature T e of the nebula -with different weights in the high and low temperature 
regions- are necessary to derive T and t 2 . In their observations, |Peimbert| ( |1967| ) found that 



the temperature derived from the Balmer discontinuity, using the intensity ratio BJ/J(H/3) 
where BJ is the value of the Balmer jump, was systematically lower than that derived from 
the [O III] 4363/5007 CELs ratio. This result has been corroborated by several subsequent 
works (see e.g. |Garcia-Kojas et al.| |2006| , and references therein) and we will refer to this 
problem as the temperature discrepancy. 



Later, [Peimbert fc Costero| (|1969|) claimed that, in the presence of temperature variations 



over the observed volume of the nebula, the gaseous abundances derived from the analysis 
of intensity ratios of CELs are underestimated if these temperature fluctuations are not 
considered. This is due to the strong dependence of the intensity of those lines on the assumed 
T e . This possibility is annoying, considering that the analysis of CELs is the standard method 
for deriving ionic abundances in ionized nebulae. On the other hand, intensity ratios of 
optical recombination lines (ORLs) are almost independent of the temperature structure of 
the nebula, and could be more appropiate to derive the "real" abundances of the nebulae. In 
the last years, several measurements of ORLs on H II regions and planetary nebulae (PNe) 
have been done, and it has been found that abundance determinations from ORLs are 
systematically higher than those obtained using CELs, independently of the ion considered 
(see e.g. |Liu| |2006| ; |Carcia-Rojas et aT| [2006| , and references therein). This is usually known 



- 3- 



as the abundance discrepancy problem. 

Thanks to the fast technological progress -related to the development of large tele- 
scopes and better detectors- and to the significant improvement in the understanding of the 
recombination theory for multielectronic atomic systems, many works have been dedicated 
to understand the origin of the abundance and temperature discrepancies. During the last 
years, it has been confirmed that both effects are real, and that they are not related neither 
to systematic biases in the observations nor in the basic atomic physics. In addition, it has 
been found that the two dichotomies can be intimately related (see [Liu et aT [2000D . 



The high surface brightness of the brightest PNe has permitted to detect many ORLs of 
heavy elements in these objects. The pioneering work of [Wyse| ( |1942| ) in NGC 7009, was the 
first in calling the attention on the different ++ abundance obtained from O II ORLs and 
[O III] CELs. Later, |Aller fc Menzel ( |I945| ) and [lbrres-Peimbert fc Peimbert| ( P77| ) also 
remarked the excess of carbon observed in some PNe when calculating the C ++ abundance 
from C II ORLs. Since the availability of the IUE satellite, carbon abundances from UV 
CELs could be determined, obtaining abundances systematically smaller than those derived 
from ORLs (see e.g. [Harrington et aL] |1980| ; |Kalei| |1986| ; |Rola fc Stasihska| |1994| ; |Peimbert| 
et |]~995|) . During the last decade many efforts have been done to understand the physical 
conditions under which the different types of spectral lines are formed in PNe. In particular, 
the use of medium-large aperture telescopes, has allowed that the catalogue of detected O II 
and C II ORLs in PNe has increased spectacularly in the last years. The works dedicated 
to the study of the abundance discrepancy and/or the temperature structure in PNe are 
fairly numerous (e.g. |Garnett fc Dinerstein] [2001| ; [Mathis et al.| |1998| ; [Liu et aE] |2001 



2006 



12004] [Feimbert et~aT| |2004] ; |Robertson-Tessi fc Garnett] |2005| ; [1'samis et aT] p004l ; |Wesson 
et al.| PUUB] , and references therein). The data collected in these works, indicate that the 
t 2 parameter shows an extraordinary range of variation in PNe, reaching very high values 
in some objects. [Liu et aL| (|2006|) , found values of t 2 ~ 0.14 and 0.18 in Hf 2-2, from the 
abundance and temperature discrepancies, respectively. On the other hand, [Rubin et al 



( |2002|) , from HST image and spectroscopy of PN NGC 7009, found that the T e map of this 
nebula was practically uniform and that t 2 throughout the nebula on the plane of the sky 
was small {t 2 a = 0.0035). In spite of the low t 2 value found, it could not completely be 
discarded the existence of temperature fluctuations throughout the line of sight that could 
account for the observed abundance discrepancies in this object (see [Liu et al] |1995| ). 



The abundance discrepancy is commonly quantified using the so-called abundance dis- 
crepancy factor (ADF), that is defined as: 



ADF(X +i ) = (X + VH + )orls/(X + 7H + )cels. 
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Liu et al. ( [200 1| ) found a correlation between the ADF and the difference between 
T e ([0 III]) and T e (Bac). Qualitatively, this correlation is consistent with the existence of 
temperature fluctuations, although these authors argue that abundance and temperature 
discrepancies are qualitatively better explained on the basis of a chemically inhomogenous or 
two-phase model. [Liu et aL| (|2000|) analized the emission line spectrum of the PN NGC 6153, 
and concluded that the nebula contains a component of the ionized gas previously unknown 
-cold and very metal-rich-. Later works have supported this model, sustaining it mainly 
in the fact that the [O III] fine structure far infrared (FIR) lines at AA 52, 88 /zm -which 
have very low excitation energies, of the order of 1000 K, and are almost insensible to the 
variations of T e - give ionic abundances comparable to the ones obtained from UV and optical 
CELs. These results are not compatible with the temperature fluctuations paradigm (see 



Liu| |2006j and references therein), but can be explained by a two-phase model (see below). 



The observational features that allow us to determine the t 2 parameter (ORLs, H I 
recombination continua) are much more difficult to be accurately measured in H II regions 
than in PNe. The low surface brightness of most of these objects has caused that, until a few 
years ago, the H II regions sample with t 2 determinations available has been limited to a few 
Galactic H II regions: the Orion nebula ( |lbrres-Peimbert et al.| |1980| ; pster brock et al.||1992 



Peimbert et aT| [Hggg: [Esteban et al] |Rubin et al.| |2PU3D , M8 ( [Peimbert et"aT| ggjgg 

Esteban et al] |1999b| ) and M17 ( Peimbert et al.||1992| ; [Esteban et "aL| |1999a| ; [Tsamis et al 



2003|) , and to a few extragalactic H II regions flGonzalez-Delgado et al.| |1994| ; [Esteban et al. 



|2002| ). In the last four years, this sample has been extended with the works of [Tsamis et al 



( |2003|) on an additional Galactic H II region (NGC 3576) and 3 extragalactic ones (1 in the 
Large Magellanic Cloud, LMC and 2 in the Small Magellanic Cloud, SMC); [Peimbert| (|2003D 
on 30 Doradus; |Peimbert et al.| ( |2005| ) on 2 NGC 6822 H II regions; and [Lopez-Sanchez et al. 



( P007| ) on the starburst galaxy NGC 5253. Aditionally, P'Dell et al] (|2003|) made a direct 
estimation of t 2 in the plane of the sky of the Orion nebula from deep WFC/HST images, and 
making some geometrical considerations found a value of £ 2 =0.028±0.006, consistent with 
previous results (t 2 between 0.020 and 0.030, see e.g. [Esteban et al.| |2004| , and references 
therein). Finally, our group has obtained and analyzed extremely deep echelle spectra of a 
sample of 8 Galactic H II regions ( Esteban et ai] 12004] ; |Garcia-Ro,jas et al.| gU0|, [2005] |2006| . 
gOUTI) . In 4 of them (M16, M20, S 311 and NGC 3603), C II and O II ORLs were detected, 
for the first time; in other two regions (Orion nebula and NGC 3576) Ne II ORLs were also 
detected for the first time. 

On the contrary, some puzzling results for the temperature fluctuations hypothesis have 
been recently found for some extragalactic H II regions. On the one hand, [Hagele et al. 



( |2006|) have measured T e (Bac) in three compact H II galaxies of the second data emission 
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release of the SDSSQ reobserved with the 4.2m William Herschel Telescope (WHT), finding 
t 2 's values close to zero for two of the objects, and a relatively large value (t 2 ~0.066±0.026) 
for the third. On the other hand, |Guseva et al.| ( |2006| , |2007| ) have made a detailed study of a 



sample of 47 low metallicity H II regions in blue compact dwarf (BCD) galaxies and and 69 
emission-line galaxies from the data release 3 of the SDSS and M101. These authors have 
measured the Balmer and Paschen jumps with the purpose of determining T e (H I). From 
Monte Carlo simulations, they fit the spectral energy distribution (SED) of the galaxies 
and find that the temperatures in the ++ zone -derived from [O III] AA(4959 5007)/A4363 
ratio- are, considering the dispersion of the data, equal to the temperature of the H + zone, 
obtained from fitting the Balmer jump and the SED. In any case, the dispersion of the data 
does not allow to discard small temperature fluctuations of about ~ 3% - 5%. 

Tsamis et a"T? ( |2003| ) proposed a chemically inhomogeneous scenario for H II regions 
similar to that proposed for PNe, but only very recently a physically reasonable model has 
been drawn ( [Tsamis fc Pequignot| |2005| ; [Stasiriska et al.| |2007|) . This scenario implies the 



presence of temperature fluctuations in the nebulae, but its resulting abundance pattern 
is not identical to the one found under the standard temperature fluctuations paradigm 
(see [Stasiriska et al.| |2007| ). Those authors postulate the presence of a gas component of 
high metallicity, and consequently low T e embedded in a less dense ambient gas, with lower 
metallicity and therefore larger temperature. According to Tsamis fc Pequignot| (|2005|) and 



IStasihska et al.| (|2007|) , this denser component could come from the resulting material of 



type II SNe that has not been mixed with the bulk of ISM and that is in pressure balance 
with the normal composition ambient gas (an scenario originally proposed by [Tenorio-Tagle 



1996|) . According to this model, these inclusions would be responsible for most the ORLs 



emission and they would not emit in CELs due to their low temperature. 

In this paper, we present a global analysis of some results of our study of a sample of 
H II regions mainly focussed on discussing the properties and behavior of the abundance 
discrepancy problem in these objects. 

The structure of this paper is as follows: in § ^ we briefly describe the data we have 
analized. In § ^ we explore the behavior of the abundance discrepancy factor (ADF) with 
respect to some nebular parameters. In § |3] we present some observational arguments against 
the two phase model in H II regions. Finally in § [5] we briefly draw the conclusions of this 
work. 



2 Sloan Digital Sky Survey 
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2. The data 

The observations were made in two runs: March 2002 and March 2003 with the Ultravi- 
olet Visual Echelle Spectrograph (UVES) at the Very Large Telescope (VLT) Kueyen unit in 
Cerro Paranal Observatory (Chile). A detailed description of the instrument and telescope 
configuration, as well as of data reduction and analysis, can be found in previous papers 
flEsteban et all [2001 ; Garcia-Rojas et al.||2004| , p005| , gOOg , g007] ; |L6pez-Sanchez et al.| p007|) . 



With the aim of presenting a homogeneous data set, we have recomputed the physical 
conditions and ionic and total abundances of some of our published results for some objects 
(NGC 3576, Orion nebula, NGC 3603) using the same atomic data and ionization correction 
factors (ICF) scheme than for the rest of the sample nebulae. Also, we have corrected some 
misprints ocurring in previous papers. In Table [l], we present the physical conditions of 
the ionized gas computed for our sample; it can be noted that only small corrections in the 
adopted T e and n e have been done with respect to the previously published ones. In Table ^ 
we present the estimated t 2 values, which are essentially the same as the previously published 
ones. Also, in Table |3] we present the compilation of the ionic abundances from ORLs. In 
Tables |] and we present the recomputed ionic and total abundances from CELs. We have 
re-calculated the electron density for the Orion nebula, NGC 3576 and NGC 3603. These 
new densities imply new values for the ionic abundances for the most density dependent 
ionic abundances i. e. + , S + , Cl ++ , Fe 3+ and Ar +3 . Changes in the adopted atomic 
data for + and S ++ (see Garcia-Rojas et aT| [2005|) have affected the determination of the 



ionic abundances and physical conditions of NGC 3576 and the Orion nebula. However, the 
combined effect of these corrections is small for NGC 3576 and NGC 3603 (less than 0.03 
dex). For the Orion nebula the effect is somewhat more important, reaching a factor of 
about 0.1 dex for + . 

Finally, we have re-computed He abundances for all the objects adopting a new set of 
ICFs: for the low ionization objects we have adopted the ICF by Peimbert fc Torres-Peimbert 



( |1977|) , whereas for the high ionization ones we have adopted the ICF by Peimbert et al. 

( TO - 



3. Behavior of the ADF with respect to some nebular parameters. 

The different hypothesis proposed to explain the abundance discrepancy predict the 
existence or absence of correlations between the ADF and some nebular properties. Although 
these correlations have been exhaustively explored for PNe, there has not been a systematic 
study for H II regions. The aim of this section is to analyze possible correlations between the 
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Table 1. Plasma diagnostics. 



Diagnostic 



Line 



M16 



M8 



M17 



M20 



1100+ 750 
iiuu_ 400 


ieoo±™° 


1200+ 1250 


ODU -220 


1050±250 


1800±800 


480+150 


240+.70 


1390±550 


1600+450 


500+.220 


320±130 


540+ > 1000 


2600±1450 


430 +>iooo 


560±390 


1370±1000 


2100±700 
2450: d 


Q7n + 630 

ziu_ 270 
>800 


ocrn + 780 
•JOU-SSO 


1120±220 


1800±350 


470±120 


270±60 



n e (cm- 3 ) [NI] (A5198)/(A5200) 
[O II] (A3726)/(A3729) 
[S II] (A6716)/(A6731) 
[Fe III] 

[CI III] (A5518)/(A5538) 
[Ar IV] (A4711)/(A4740) 
n c (adopted) 



To (K) 


[N II] (A6548+A6583)/(A5755) a 


8450+270 


8470±180 


8950±380 


8500±240 




[S II] (A6716+A6731)/(A4069+A4076) 


7520±430 


7220±300 


7100±750 


6950 ±350 




[0 II] (A3726+A3729)/(A7320+A7330) a 


8260+.400 


8700±350 


8750±550 


8275±400 




T c (Low) 


8350+.200 


8500±150 


8870±300 


8400±200 




[O III] (A4959+A5007)/(A4363) 


7650+.250 


8090±140 


8020±170 


7800±300 




[Ar III] (A7136+A7751)/(A5192) 




7550±420 


8380±570 


8730±920 b 




[S III] (A9069+A9532)/(A6312) 


8430±450 


8600±300 c 


8110+400 


8300±400 




To (High) 


7850±220 


8150+120 


8050±150 


7980±250 




T c (O II+III) 


8180±300 


8570±200 


8200±200 


8230±350 




He II 


7300±350 


7500±200 


7450±200 


7650±300 




Balmcr discontinuity 


5450±820 


7100±1100 




6000±1300 




Paschen discontinuity 


7200±1300 


7750+900 


6500+1000 


5700+1300 






NGC 3576 


ORION 


NGC 3603 


S 311 



n c (cm" 3 ) [N I] (A5198)/(A5200) 
[O II] (A3726)/(A3729) 
[S II] (A6716)/(A6731) 
[Fe III] 

[CI III] (A5518)/(A5538) 
[Ar IV] (A4711)/(A4740) 
n c (adopted) 



1900±400 

1300+ 500 

3100± 1300 

35001?°° 

4500+.?°°° 

2300±300 



1700±600 

6300+ 2800 

6500+ 2000 
oouu_ 1200 

9300+2700 

94001^°° 

68001J™ 

7800±600 



4000: d 
3000+1000 

4150 +3350 

1330: d 
5600+ 3900 

ODUU_ 2400 

< 4000 
3400±850 



-,+260 
-200 



590Z 
260+11 ! ; 
360 4 



^140 
-120 

390±220 
550+ 650 

310±80 



To (K) 


[N II] (A6548+A6583)/(A5755) a 


8500±200 


10150±350 


11050±800 


9500±250 




[S II] (A6716+A6731)/(A4069+A4076) 


8400l 3 go 


9050±800 


iiuou_ 2050 


7200+ 750 




[O II] (A3726+A3729)/(A7320+A7330) a 


8050±450 


8700±500 


12350±1250 


9800±600 




T e (Low) 


8400±200 


9600±300 


11400±700 


9550±250 




[O III] (A4959+A5007)/(A4363) 


8500±50 


8300±40 


9060±200 


9000±200 




[Ar III] (A7136+A7751)/(A5192) 


8600l«° 


8300±400 




88001™ 




[S III] (A9069+A9532)/(A6312) 


8750±550 


9700+ 800 


8800 ±500 c 


9300±350 c 




To (High) 


8500±50 


8320±40 


9030±200 


9050±200 




T (O II+III) 


8500±50 


8430±130 


9600±200 


9600±450 




He II 


6800±400 


7950±200 


8480±200 


8750±500 




Balmcr discontinuity 


6650±750 


7900±600 




9500±900 




Paschen discontinuity 


6700+900 


8100±1400 


6900±1100 


8700±1100 


a The recombination contribution to the auroral lines has been taken into account. 



b The [Ar in] A7751 line is severely blended with a telluric line. 
C [S III] A9069 or A9530 affected by atmospheric absorption bands. 

d Colons indicate very high uncertainties. These values has not been taken into account in the adopted average. 
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Table 2. t 2 parameter. 



Method 


M16 


M8 


M17 


M20 


0+ (R/C) 




0.031±0.017 


0.109: 


0.032±0.020 


0++ (R/C) 


0.046±0.007 


0.045±0.005 


0.034±0.005 


0.038±0.016 


C++ (R/C) 




0.035±0.005 






Ne++ (R/C) 










He+ 


0.017±0.013 


0.046±0.009 


0.027±0.014 


0.017±0.010 


Bac/Pac-LEC 


0.045±0.014 


0.022±0.015 


0.035±0.021 


0.049±0.019 


Adopted 


0.039±0.006 


0.040±0.004 


0.033±0.005 


0.029±0.007 


Method 


NGC 3576 


ORION 


NGC 3603 


S 311 


0+ (R/C) 




0.063±0.029 






0++ (R/C) 


0.038±0.006 


0.020±0.002 


0.042±0.009 


0.040±0.008 


C++ (R/C) 




0.039±0.011 






Ne++ (R/C) 


036+ 014 

U.UOD_ Q24 


0.034±0.014 






He+ 


0.023±0.019 


0.022±0.002 


0.032±0.014 


0.034±0.010 


Bac/Pac-LEC 


0.037±0.012 


0.016±0.014 


0.056±0.023 


0.009: 


Adopted 


0.038±0.009 


0.022±0.002 


0.040±0.008 


0.038±0.007 



Table 3. Heavy element ionic abundances from ORLs. 



12 + log(X+7H+) 
Ion M16 M8 M17 M20 



0+ ••• 8.53±0.06 •■• 8.62±0.07 

0++ 8.30±0.04 8.23±0.02 8.68±0.02 8.00±0.18 

C++ 8.40±0.03 8.30±0.02 8.68±0.03 8.18±0.05 
Ne++ 

Ion NGC 3576 ORION NGC 3603 S 311 



o+ 

0++ 8.62±0.05 8.57±0.01 8.71±0.05 8.08±0.03 

C++ 8.45±0.06 8.34±0.02 8.48±0.07 8.00±0.04 

Ne++ 7.90±0.14 7.95±0.09 
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Table 4. Heavy element ionic abundances from CELs. 



12 + log(X+VH+) 



Ion 



M16 



M8 



M17 



M20 



0.000 



0.039 



0.000 



0.040 



0.000 



0.033 



0.000 



0.029 



N° 


6 


15±0 


06 


6 


33±0 


07 


5 


81±0 


05 


5 


99±0 


06 


5.57±0.07 


5.70±0.07 


5 


90±0 


07 


6 


03±0.08 


N+ 


7 


71±0 


05 


7 


88±0 


06 


7 


50±0 


03 


7 


67±0 


04 


6.82±0.10 


6.94±0.10 


7 


55±0 


04 


7 


67±0.05 


o° 


7 


23±0 


05 


7 


40±0 


06 


6 


40±0 


03 


6 


57±0 


04 


6.87±0.07 


7.00±0.07 


6 


60±0 


05 


6 


72±0.06 


o+ 


8 


47±0 


08 


8 


66±0 


09 


8 


39±0 


06 


8 


58±0 


07 


7.84±0.09 


7.98±0.09 


8 


46±0 


07 


8 


59±0.08 


o++ 


7 


85±0 


07 


8 


18±0 


10 


7 


86±0 


03 


8 


18±0 


07 


8.41±0.04 


8.67±0.06 


7 


67±0 


08 


7 


90±0.10 


Ne++ 


7 


01±0 


07 


7 


38±0 


10 


6 


95±0 


05 


7 


30±0 


07 


7.64±0.04 


7.93±0.07 


6 


55±0 


09 


6 


80±0.11 


S+ 


6 


32±0 


05 


6 


49±0 


06 


5 


93±0 


04 


6 


10±0 


07 


5.44±0.05 


5.56±0.06 


6 


17±0 


05 


6 


29±0.06 


S++ 


6 


84±0 


06 


7 


22±0 


10 


6 


89±0 


03 


7 


25±0 


07 


6.90±0.04 


7.19±0.06 


6 


79±0 


06 


7 


09±0.10 


C1+ 


1 


77±0 


05 


4 


91±0 


07 


4 


53±0 


04 


4 


66±0 


06 


3 Q5+ 09 

o,aJ -0.12 


4 06+ 09 


i 


75±0 


05 


4 


85±0.07 


C1++ 


5 


04±0 


06 


5 


36±0 


08 


5 


02±0 


04 


5 


32±0 


06 


5.04±0.04 


5.29±0.06 


1 


99±0 


07 


5 


21±0.08 


Cl 3 + 


























3.15: 


3.35: 












Ar++ 


6 


25±0 


05 


6 


53±0 


08 


6 


21±0 


03 


6 


48±0 


05 


6.35±0.04 


6.57±0.06 


6 


17±0 


06 


6 


36±0.08 


Ar 3+ 


3 


89±0 


22 


4 


23±0 


23 


3 


69±0 


09 


4 


01±0 


10 


4 15+ ' 12 


4 ' 4Z -0.18 


4 


01±0 


18 


4 


24±0.19 


Fe+ 




4.62: 






4.78: 






4.61: 






4.77: 




4.05: 


4.17: 




4.51: 






4.62: 


Fe++ 


5 


07±0 


04 


5 


41±0 


08 


5 


58±0.04 


5 


91±0 


06 


5.24±0.06 


5.51±0.08 


5 


23±0 


10 


5 


47±0.12 


Fe 3 + 









































12 + log(X+VH+) 



NGC 3576 



ORION 



NGC 3603 



S 311 



Ion 


t 2 = 0.000 


t 2 = 0.038 


t 2 = 0.000 


t 2 = 0.022 


e 


= 0.000 


t 2 


= 0.040 


t 2 = 0.000 


t 2 = 0.038 


N° 








5.69±0.06 


5.76±0.06 


5 


55±0 


11 


5 


65±0 


11 


5.74 ± 0.06 


5.88 ± 0.07 


N+ 


7.09±0 


06 


7.25±0.07 


6.97±0.09 


7.03±0.09 


6 


45±0 


07 


6 


55±0 


07 


7.25 ± 0.05 


7.38 ± 0.06 


0° 


6.35±0 


04 


6.51±0.06 


6.22±0.05 


6.28±0.05 


6 


32±0 


09 


6 


42±0 


09 


6.74 ± 0.06 


6.87 ± 0.06 


o+ 


8.21±0 


07 


8.38±0.08 


7.93±0.15 


8.00±0.15 


7 


42±0 


11 


7 


52±0 


11 


8.26 ± 0.07 


8.40 ± 0.08 


o++ 


8.35±0 


03 


8.63±0.08 


8.42±0.01 


8.58±0.03 


8 


42±0 


05 


8 


68±0 


08 


7.81 ± 0.04 


8.05 ± 0.06 


Ne++ 


7.61±0 


09 


7.91±0.12 


7.67±0.07 


7.84±0.07 


7 


72±0 


05 


8 


00±0 


08 


6.81 ± 0.05 


7.07 ± 0.07 


S+ 


5.75±0 


08 


5.91±0.09 


5.44±0.06 


5.50±0.06 


5 


08±0 


10 


5 


17±0 


10 


6.03 ± 0.05 


6.15 ± 0.06 


S++ 


6.86±0 


10 


7.17±0.10 


6.95±0.04 


7.12±0.05 


6 


83±0 


01 


7 


11±0 


09 


6.68 ± 0.07 


6.95 ± 0.09 


C1+ 


4.13±0 


08 


4.26±0.08 


3.80±0.11 a 


3.85±0.11 a 


3 


46±0 


07 


3 


54±0 


07 


4.56 ± 0.05 


4.67 ± 0.05 


C1++ 


4.95±0 


06 


c 9-1+0.10 
°- zi -0.08 


5.13±0.02 


5.28±0.02 


5 


07±0 


05 


5 


31±0 


08 


4.85 ± 0.05 


5.08 ± 0.05 


Cl 3 + 


3.21±0 


07 


3.42±0.09 


3.66±0.12 


3.78±0.12 


3 


86±0 


04 


4 


06±0 


07 






Ar++ 


6.34±0 


05 


6.57±0.08 


6.37±0.05 


6.50±0.05 


6 


35±0 


04 


6 


56±0 


07 


6.08 ± 0.04 


6.28 ± 0.06 


Ar 3+ 


4.20±0 


07 


4.48±0.10 


4.59±0.03 


4.75±0.04 


4 


85±0 


06 


5 


11±0 


08 


, 4 2+0.18 

-0.30 


3 fifi+O-18 
J.DD_0 .30 


Fe+ 


4.61: 




4.76: 


4.54: 


4.60: 




4.04: 






4.13: 




4.25: 


4.37: 


Fe++ 


5.57±0.05 


5.85±0.09 


5.43±0.06 


5.59±0.06 


5 


24±0 


06 


5 


50±0 


09 


5.05 ± 0.06 


5.30 ± 0.08 


Fe 3+ 


°- "--0.29 


°- MO -0.29 


5 65+ 019 
ODO -0.30 


r- 70+O.I9 
°- ,<5 -0.30 



















a These values were misprinted in Esteban et al. (2004) 
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Table 5. Total abundances. 



12 + log(X+VH+) 



M16 



M8 



M17 



M20 





t 2 = 0.000 


t 2 = 0.039 


t 2 = 0.000 


t 2 = 0.040 


t 2 = 0.000 


t 2 = 0.033 


t 2 = 0.000 


t 2 = 0.029 


He 


11.01±0.02 


10.97±0.02 


10.87±0.01 


10.85±0.01 


10.97±0.01 


10.97±0.01 


10.95±0.06 


10.92±0.06 


C a 


8.85±0.10 


8.85±0.10 


8.61/8.69±0.09 


8.69/8.69±0.09 


8.77±0.04 


8.77±0.04 


8.66±0.11 


8.66±0.11 


N 


7.84±0.06 


8.07±0.12 


7.72±0.03 


7.96±0.06 


7.62±0.12 


7.87±0.13 


7.67±0.05 


7.83±0.07 





8.56±0.07 


8.78±0.07 


8.51±0.05 


8.73±0.05 


8.51±0.04 


8.75±0.05 


8.53±0.06 


8.67±0.07 


o b 


8.70±0.06 


8.70±0.06 


8.71±0.05 


8.71±0.05 


8.89±0.04 


8.89±0.04 


8.71±0.09 


8.71±0.09 


o c 


8.81±0.07 


8.81±0.07 


8.71±0.04 


8.71±0.04 


8.76±0.04 


8.76±0.04 


8.69±0.10 


8.69±0.10 


Ne d 


7.86±0.15 


8.08±0.17 


7.81±0.12 


8.03±0.13 


7.74±0.07 


8.01±0.09 


7.83±0.16 


7.97±0.18 


S 


6.96±0.05 


7.29±0.08 


6.94±0.03 


7.28±0.06 


7.01±0.04 


7.33±0.06 


6.88±0.05 


7.12±0.09 


Cl e 


5.23±0.04 


5.49±0.07 


5.14±0.04 


5.41±0.06 


5.08±0.04 


5.32±0.06 


5.19±0.05 


5.37±0.06 


Ar 


6.70±0.06 


6.84±0.08 


6.52±0.05 


6.69±0.06 


6.39±0.06 


6.59±0.07 


6.65±0.07 


6.70±0.09 


Fe f 


5.17±0.11 


5.53±0.13 


5.69±0.09 


6.04±0.10 


5.87±0.13 


6.22±0.15 


5.31±0.13 


5.56±0.15 


Fes 


5.20±0.06 


5.51±0.07 


5.62±0.04 


5.94±0.06 


5.27±0.06 


5.52±0.08 


5.31±0.09 


5.52±0.10 


12 + log(X+7H+) 



NGC 3576 



ORION 



NGC 3603 



S 311 





t 2 = 0.000 


t 2 = 0.038 


t 2 = 0.000 


t 2 = 0.022 


t 2 = 0.000 


t 2 = 0.040 


t 2 = 0.000 


t 2 = 0.038 


He 


10.97±0.03 


10.96±0.04 


10.95±0.004 


10.95±0.004 


10.99±0.01 


10.99±0.01 


10.99±0.02 


10.97±0.02 


C a 


8.61±0.08 


8.61±0.08 


8.43±0.04 


8.43±0.04 


8.51±0.07 


8.51±0.07 


8.37±0.10 


8.37±0.10 


N 


7.57±0.06 


7.81±0.07 


7.73±0.09 


7.87±0.09 


7.62±0.13 


7.89±0.14 


7.43±0.06 


7.61±0.07 


O 


8.59±0.03 


8.82±0.07 


8.54±0.03 


8.68±0.04 


8.46±0.05 


8.71±0.07 


8.39±0.05 


8.56±0.06 


O b 


8.74±0.06 


8.74±0.06 


8.71±0.03 


8.71±0.03 










o c 


8.82±0.06 


8.82±0.06 


8.67±0.03 


8.67±0.03 


8.72±0.05 


8.72±0.05 


8.57±0.05 


8.57±0.05 


Ne d 


7.85±0.10 


8.11±0.12 


7.79±0.07 


7.94±0.07 


7.76±0.08 


8.03±0.11 


7.79±0.13 


7.98±0.14 


S 


6.92±0.10 


7.23±0.10 


7.04±0.04 


7.23±0.04 


7.04±0.05 


7.36±0.08 


6.77±0.06 


7.02±0.08 


CP 


5.02±0.05 


5.27±0.07 


5.14±0.04 


5.29±0.04 


5.11±0.05 


5.34±0.07 


5.03±0.06 


5.22±0.07 


Ar 


6.44±0.06 


6.65±0.09 


6.50±0.06 


6.63±0.06 


6.37±0.06 


6.58±0.08 


6.43±0.06 


6.56±0.07 


Fe f 


5.89±0.21 


6.23±0.10 


5.96±0.16 


6.18±0.14 


6.14±0.16 


6.53±0.19 


5.17±0.11 


5.44±0.13 


FeS 


5 07+011 

°' M '-0.15 


6.22±0.12 


5.88±0.12 


6.02±0.12 


5.27±0.06 


5.52±0.09 


5.11±0.06 


5.35±0.08 



a For M8 and Orion nebula, ICFs from C II] A2525 line/ICF from photoionization models by Garnctt et al. (199E) 
b O+/H+ and O++/H+ from O I and O n ORLs. 
c O+/H+ from [O n] CELs and i 2 >0.00. 0++/H+ from O n ORLs. 
d ICF from 



Peimbert & Costerc 



(1969) for M17, NGC 3576 and NGC 3603. ICF from Torres- Peimbert & Peimbert (1977) for 



M16, M8, M20 and S 311. 

C C1/H=C1+/H+ + C1++/H+ for M16, M8, M20, and S 311. C1/H=C1+/H 4 
Orion nebula and NGC 3603. 



C1++/H+ + Cl+ 3 /H+ for M17, NGC 3576, 



'ICF from 



Rodriguez & Rubin (2005) equation (2) 



SFe/H=Fe+/H+ + Fe++/H+ for M16, M8, M17, M20, NGC 3603 and S 311. Fe/H=Fe+/H+ + Fe++/H+ + Fe+ 3 /H+ for 
NGC 3576 and Orion nebula. 
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ADF and different nebular properties in H II regions (from our data and others available in 
the literature), with the purpose of discarding possible systematic errors and verifying the 
consistency of some hypothesis proposed to explain the abundance discrepancy. 



Hereafter, the ADF found by [1 samis et aT| (|2003|) for the SMC N11B H II region has 



been omitted in most of the figures and the discussion, because the ADF value for this region 
is much larger than for the rest of the H II regions and may be considerably overestimated. 
For this object, [l'samis et al.| ( |2003| ) made an attempt to correct the intensity of the RLs 



of multiplet 1 of O II for the presence of stellar absorption features, mainly caused by B 
stars -that have a strong absorption O II spectra- in the area covered by the slit. In any 
case, this effect, which can be important in extragalactic objects, can only be corrected 
in the appropriate form if the spectral features of the stars are resolved, or if synthetic 
spectra are available, and this is not the case in this particular object. In fact, it is not 
possible to distinguish between the nebular emission and stellar absorption features in the 
spectra of [Tsamis et al.| ( |2003|) of N11B. It is important to emphasize that this effect must be 



investigated whenever is possible to perform a suitable analysis when computing abundances 
from ORLs in extragalactic H II regions. 

For all the fits, we have computed their Spearman's correlation coefficiens. These coeffi- 
cients are the most appropriate ones for non gaussian distributions, as well as the two-sided 
significance of its deviation from zero; the significance is a value in the interval [0.0, 1.0], 
and a small value indicates a significant correlation. 

One important result of our data compilation for H II regions is that the ADF is fairly 
constant for these objects and of the order of 2. This is a crucial difference with the behavior 



shown by this factor in PNe, where it can vary from values of 1 up to 20 (see § L2 for a 
further discussion and references). 

In Figure p] to |3] we have represented the ADF with respect to different nebular pa- 
rameters, and in Table || we include the slope of the fits, their corresponding Spearman's 
coefficients, r, and the significance of the correlation, P, for each of the aforementioned dia- 
grams. As it can be seen in Table [] the correlation coefficients for all the fits are very low and 
they have quite poor statistical significance. As in the rest of the cases, Figure [3] shows that 
the ADF seems to be not related to ++ /H + nor + /0 ++ , at least within the observational 
uncertainties. The metallicity dependency of the t 2 parameter was suggested by |Garnett 



( |1992|) , who found that photoionization models could produce temperature fluctuations sim- 
ilar to the observed ones for nebulae with T c < 9000 K, increasing for colder (more metallic) 
nebulae. However, the observed t 2 does not seem to depend on metallicity ,both for H II 
regions and PNe. For PNe, several authors have found that the large range of abundance 
discrepancies found might be mainly due to the fact that ADF increases monotonically with 
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Fig. 1. — ADF(0 ++ ) vs. the degree of ionization (upper panel) and ++ /H + determined 
from CELs (lower panel). The Galactic H II regions of our sample are represented by filled 
dots; the rest of symbols represents data from the literature for extragalactic H II regions: 
filled triangles: |L6pez-Sanchez et al.| ( |2007| ); filled squares: |Tsamis et al.| ( |20Q3| ) ; empty 
Fsteban et al. 



dots: 



; empty triangles: [Peimbert fc Peimbert] (|2005|); empty squares: 



Peimbert] tfM)§. 
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metallicity (see e.g. |Liu et aL||200U| , |2U04| ). In the case of H II regions, as it can be observed 
in Figure the ADF seems to remain relatively constant with the abundance, determined 
either from CELs or ORLs, and within the uncertainties. We have also represented the ADF 
vs. the ratio between high and low ionization temperatures, finding that there is, again, no 
correlation within the observational errors (see Figure § and Table §); this indicates that 
large scale variations of T e due to natural temperature gradients throughout a nebula, do 
not seem to be related to the abundance discrepancy. As it can also been seen in Figure [| 
we have verified that the ADF does not seem to depend on the assumed temperature, which 
discards systematic effects in the determination of abundances from CELs, at least at the 
precision of our data. 



3.1. Testing the effect of gas kinematics. 

Several authors have proposed that mechanical energy deposited by shocks is a mech- 
anism to increase the heating of nebulae and a possible source of temperature fluctuations 
(e.g. [Peimbert| |1995|; [Peimbert et al.||i991|). One of the manifestations of the shocks consists 



of line broadening or line splitting due to the presence of different kinematical components. 
Hence, a way to investigate whether the ADF can be related to excitation by shocks is rep- 
resenting it with respect to the full width at half maximum (FWHM) of different emission 
lines. In Figure | we represent the ADF vs. the FWHM of H/3, [O III] A4959 and [N II] A6548. 
In order to extend the FWHM baseline, we have introduced the points corresponding to 3 
slit positions of NGC 5253 analyzed by |L6pez- Sanchez et al.| ( f2007|) . We have corrected the 



FWHM by the instrumental width. It is possible to note that FWHM are much larger for the 
extragalactic objects due to large scale gas movements, that are not observed in the Galactic 
H II regions. With the available data, we cannot find a clear relation between the ADF and 
the represented FWHMs. It is obvious that the paucity of points with large FWHM does 
not permit to conclude any strong statement about this possible relation. Therefore, it is 
necessary to increase the number of extragalactic objects with good determinations of both 
quantities for making a suitable exploration. 



3.2. Testing the effect of dust grain heating. 

Another possible effect that can affect the ADF was proposed by [Stasihska fc Szczerba 
( |2001|) and consists of photoelectric heating of dust grains in the zones near the ionizing star 
of a PN that increase T e in the inner zones producing large t 2 values. This is the case if 
the nebula is inhomogeneous in density, because the fluctuation in the ionization parameter 
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Fig. 2. — ADF(0 ++ ) vs. total abundance derived from CELs (upper panel) and ORLs 
(lower panel). Symbols are the same that in Figure |I[ 
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Fig. 3. — ADF(0 ++ ) vs. T e (low) /T e (high) (upper panel) and vs. T e (high) (lower panel). 
Symbols are the same that in Figure ffl. 
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Fig. 4.— ADF(0 ++ ) vs. FWHM of H/3 (upper panel), [0 ill] A4959 (medium panel) and 
[N il] A6548 (lower panel). Symbols are the same that in Figure |1[ 
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induced by the variation of density, produce t 2 due to differential dust and hydrogen heating. 
In the case of an homogeneous density, dust produces a temperature gradient but not real 
temperature fluctuations flStasiriska fc Szczerbaj [200 1|) . A way to verify if the effect of dust 
grains is significant in the observed ADF is to look for a correlation between the effective 
temperature, T e ff, of the ionizing star of the nebula -that gives an indication of the hardness 
of the radiation field- and the ADF found in the associated nebula. [Robertson- Tessi 
Garnett| ( |2005|) obtained that for PNe, the ADF was not correlated with T e ff of the ionizing 
star. For H II regions this is a complicated task because they are generally ionized by an 
OB association. In Table [7| we include the T e ff of the main ionizing star (the hottest one) of 
each studied H II region, and it is clear that such quantity seems to be not correlated with 
the ADF of the associated nebula. 



3.3. The ADF and recombination continuum/CELs temperature difference. 



The standard temperature fluctuations hypothesis predicts that the ADF is related to 
the difference between T e (CELs) and T e (H I) ( [Feimbert fc (Jostero| |1969| ; [forres-Peimbert 



et al.| |i~980[) . Indeed, [Liu et ak| (|2000|) found a strong correlation between these two dis- 
crepancies from data for PNe that cannot be reproduced by the temperature fluctuation 
paradigm. The results of our work seem to indicate that there is also a slight correlation for 
H II regions, but with a flatter slope than for PNe (see Figure |5|). In Figure |5] we have also 
included the results obtained for 2 extragalactic H II regions: 30 Doradus ([Peimbert| |2003| ) 
and NGC 23630 flGonzalez-Delgado et all |199l |Esteban et aT| |2002|) . In the case of H II 
regions, contrary to what happens in PNe, the relation between ADF and the temperature 
differences seem to be consistent with the standard temperature fluctuations hypothesis, and 
they are translated into moderate and similar values of the t 2 parameter, without reaching 
the extreme values found in PNe. The fit obtained for our data is: 



log[ADF(0 ++ )] = (0.184 ± 0.022) + (0.609 ±0.208) (2) 

xl0- 4 [T e ([OIII])-T e (HI)], (3) 

with a correlation coefficient of r= 0.575. 



3 For this object the ADF has been obtained from the analysis by Esteban et al. ( 2002 ), and temperatures 
were obtained from Gonzalcz-Delgado et al. (1994). We have included these data because they correspond 
to the same slit position, covering a very similar volume of nebula. 
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Table 6. Parameters of the linear least-squares fits of the data represented in Figures [I] 

to|. 



ADF(0 ++ ) vs. 


Slope 


r 


P 


0+/0++ 


0.099 


0.204 


0.411 


12+log(0 ++ )/H+ 


-0.164 


-0.198 


0.420 


12+log(0/H) C£L 


-0.049 


0.064 


0.768 


12+log(0/H) Ofl L 


0.088 


0.283 


0.283 


T c (Low)/T e (High) 


-0.379 


0.225 


0.378 


T e (High) 


2.21xl0~ 6 


-0.113 


0.646 



Table 7. Comparison of the ADF with the spectral type of the main ionizing sources of 

each of our sample H II regions. 



Object a 


Star/Cluster 


Spectral type b 


ADF(0++) 


NGC 3603 


NGC 3603/HST-38/40/A2/16 


03V 


1.9 




NGC 3603/HST-42/A3 


02-3111 




M16 


HD168075 


04V 


2.8 


M17 


Kleinmann's star 


04V 


2.1 


M20 


HD164492 


06V 


2.1 


Orion 


^Ori C 


06peV 


1.4 


S 311 


HD 64315 


06eV 


1.9 


M8 


H36 


07.5V 


2.0 


NGC 3576 


Obscured cluster 


OB association 


1.8 



a Ordered from earlier to later spectral type. 

b From the Galactic O star catalog by [Maiz-Apellaniz et al.| ( |2004| ) . 
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£000 4000 6000 

T e ([0 III])-T e (H I) (K) 



Fig. 5. — ADF(0 ++ ) vs. the difference between T e (CELs) and T e (H i). Solid line corresponds 
to the fit found by [Liu et al.] (|2000|) for PNe, whereas pointed line corresponds to the fit to our 
data of H II regions. We have added the data of two extragalactic H II regions: 30 Doradus 
( |Feimbert| p003|) and NGC 2363 flGonzalez-Delgado et al.] |T99l ; [Fsteban et g]gQDg ). 
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Finally, as we previously pointed out (see |Garcia-Kojas et aT. 2007| ), we do not find the 
correlation between the ADF and the electron density, n c , obtained by [Robertson-'l'essi "Sz 
Garnett) (|2005|) for a PNe sample. 



4. Objections to the two phase model. 

The fact that H II regions as well as PNe are photoionized nebulae does not imply that 
photoionization is responsible for all the observed properties of these objects. A common 
argument against temperature fluctuations in H II regions is based on the difficulties of pho- 
toionization models to reproduce t 2 values of the order of those estimated for most objects. 
However, photoionization models might not be sufficiently realistic, and/or additional en- 
ergy sources could be necessary to explain the discrepancies between observations and model 
predictions (see |Viegasj p002j , and references therein). Several authors have proposed the ex- 



istence of a plasma component -colder and metal-rich- (the so-called two-phase model) to 
explain the dichotomy between abundance and temperature determinations in PNe, where 
the standard temperature and/or density fluctuations hypothesis do not seem to explain the 
available observations (|Liu et al.| |2001| ; [Tsamis et al.|p004l ). As we have already commented 



in § U, |Tsamis fc Pequignot| ( |2005|) have proposed a similar explanation for the abundance 
discrepancy in H II regions although there are not direct observational evidences of the 
colder and metal-rich inclusions they propose. In the following, we are going to present some 
objections to this model, based on the results obtained in this work. 

In Table ^] we present values of the ADF computed for different ions in Galactic and 
extragalactic H II regions. The table includes all the data available up to date, including 
our own results and others from the literature. As we commented before, the ADF for ++ 
is quite similar in all the objects -Galactic or extragalactic- and also similar to the ADFs 
for the other ions: C ++ , + , and Ne ++ . Comparing the results for H II regions of different 
galaxies, we find that the ADF seems to be independent of the morphological type, mass, 
metallicity or even the star formation history of the host galaxy (see also [Lopez-Sanchez 



|et al.|p007j ). In the case of the Galactic objects, the similarity of the ADFs suggests that the 
process that produces this phenomenon is independent of the conditions and properties of 
the Galactic disk, at least along the range of Galactocentric distances covered by our sample. 
Lopez-Sanchez et al.| Q2007] ) have discussed this fact in a paper that presents the ADF values 



3 It should be remarked that the hypothesis involved in two phase scenario are in the classical temperature 
fluctuations paradigm are not of the same order, since the two phase scenario actually gives a natural 
explanation for the origin of the temperature fluctuations. 
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for several zones of NGC 5253, a dwarf Wolf-Rayet galaxy. 



One of the observational arguments against pure temperature fluctuations in PNe by 
several authors is the agreement between total abundances derived from optical and FIR 
CELs in PNe (see e.g. [Liu et aT| [200 1| ; [Tsamis et al.|[2004|) . Due to their low excitation 
energy, E ex ~ 10 3 K, the fine structure infrared line emissivities have a very weak dependency 
with T e , therefore they must be insensible to the uncertainties introduced by the presence 
of temperature fluctuations and, in the case of the existence of such fluctuations, they must 
give values of abundance similar to that obtained from ORLs. In the case of PNe, this 
comparison is usually feasible due to the small angular size of these objects, both in the 
optical range (using scanning techniques to cover the whole nebula, see e.g. [Liu et al] p000| ), 
and in FIR and UV (see e.g. [Liu et at] |2000| , |2001|) . However, we think that this comparison 
is not so feasible in the case of extended Galactic H II regions. These objects cover much 
larger sky areas than PNe, so large uncertainties may arise due to ionization stratification 
when comparing ionic abundances derived from narrow slit optical spectroscopy and FIR 
observations obtained with space-borne facilities, that use larger apertures. In a previous 
paper ( Garcia- Roj as et~aT. 2006 ) we have showed the difficulties found when comparing our 
derived total abundances for three objects of our sample (M17, NGC 3603 and the Orion 
nebula) with those derived from FIR observations ( [Simpson et al.||l995| ). 



4.1. Dependence of the ADF on the excitation energy and the critical density. 

By definition, in the temperature fluctuations formulation, the abundance discrepancy 
must be related to the excitation energy, E ex of the upper level in which the line originates 
(see basic equations in |Peimbert] |1967|) ; whereas in the presence of density fluctuations the 



abundance discrepancy is maximized if the abundance has been calculated from a CEL 
coming from a level with low critical density n c . This last case implies that the ADF and n c 
are inversely related. [Liu et af] ( |2000| , [200 1|) , comparing different abundance determinations 



in PNe (UV, optical and FIR CELs), showed that the ADF was not related to E ex nor 
to n c . We cannot make comparisons between different types of lines (UV, optical or FIR) 
of the same ion, since we do not have observations of the same zones in other spectral 
range than the optical one; however, we can compare the ADFs obtained for Galactic and 
extragalactic H II regions, with E ex and n c of the higher level of the main CELs of each ion. 
In the case of C ++ the comparison between UV and optical data could be problematic due 
to the different volumes of nebula covered by the slits. For NGC 5253, we are sure that 
optical and UV observations cover almost the same zones, as was pointed out by [Lopez- 



|Sanchez et al.| ( |2007|) . In the case of M8 and the Orion nebula, the observations do not cover 
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Table 8. Observed ADFs in Galactic and extragalactic H II regions. 



ID 


ra e (cm 3 ) 


T e ([0 III]) (K) 


ADF(0++) ADF(Q+) 


ADF(C++) a 


ADF(Ne++) 


Ref. b 


Galactic H II regions 


M16 


1120 


7650 


2.8 






1 


M8 


1800 


8150 


2.0 1.4 


1.8 




1 




1750 


8050 


2.0 


2.1 




2 


M17 


470 


8050 


2.1 4.2: c 






1 




860, 520 


8120, 8210 


1.8, 2.2 






3 




600-1500 


8200 


2.1 






4 


M20 


270 


7980 


2.1 1.5 






1 


NGC 3576 


2300 


8500 


1.8 




2.0 


1 




1300-2700 


8850 


1.8 






4 


Orion 


7800 


8320 


1.4 1.6 C 


1.9 


1.9 


1 




4000, 5700 


8300, 8350 


1.3, 1.5 


2.5, 2.2 




5 


NGC 3603 


3400 


9030 


1.9 






1 


S 311 


310 


9050 


1.9 






1 


Extragalactic H II regions 


LMC 30 Dor 


370-1800 


10100 


2.0-2.7 






4 




316 


9950 


1.6 






6 


LMC N11B 


80-1700 


9400 


4.9-8.2 






4 


SMC N66 


50-3700 


12400 


2.3 






4 


NGC 604 


< 100 


8150 


1.6 






7 


NGC 2363 


360-1200 


15700 


2.2 






7 


NGC 5461 


300 


8600 


1.9 






7 


NGC 5471 


220-1150 


14100 


1.6 






7 


NGC 6822 V 


175 


11900 


1.9 






8 


NGC 5253 A 


580 


12100 


1.8 


2.6 




9 


NGC 5253 B 


610 


12030 


1.5 


2.1 




9 


NGC 5253 C 


370 


10810 


1.9 


2.1 




9 



a CELs abundances for C^ 



et al 



1992 



Kobulnicky et al 



were obtained from IUE data of C in 
for M8 



1997 



AA1906+1909 line. (Peimbcrt et al 
the Orion nebula and NGC 5253, respectively). 



1993b: Walter 



b l) Th i s work ; 2) 
Peimbert] ( ^003[ ); 7) 

c Dubious value. 



jEstcban et alj jl999 
pistcban ct^ff (2002 



|l999b| ); 3 ) |Esteban ct ah (|1999a 

pOO^); 8) Peimbcrt ct al. (200f); 9) [Lopez-Sanchez ct alj fcOOjJ 



4) Tsamis ct al. j iooij ) ; 5) 



Esteban et al 



(199J); 6) 
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the same zone, but we have compared the ionization correction factor for C ++ obtained 
from computing C + /H + ratio from IUE observations of the UV C n] A2326 line ( |Peimbert 



et al.| 1993b ; [Walter et al.| 1992j , for M8 and the Orion nebula, respectively) and by using 
the ionization correction factor (ICF) of |Garnett et al.| ( p. 999] ). From these data we have 
obtained ICFs of 2.48 and 2.45 for M8, and 1.20 and 1.24 for the Orion nebula, suggesting 
that there are not significant changes in the ionization structure of the observed zones in the 
optical and UV. In Figure |6] it can be noted that, taking into account the uncertainties, a 
slight correlation (r=0.7, P=0.188) between the ADF and E ex is apparent. There is also an 
apparent correlation (r—0.7, P=0.188) with n c , the opposite to that predicted by density 
fluctuations, which suggests that density fluctuations should play a much less important role 
than temperature fluctuations in H II regions. These results provide an additional argument 
that the mechanism that produces the abundance discrepancy in H II regions should be 
different to that proposed for PNe and, therefore, we can not make the same considerations 
when working with PNe and H II regions. In Table ^|, we show the values of ADF, E ex and 
n c adopted for each ion. The values we have adopted are the weighted mean of the ADFs 
computed for each ion in different objects. We distinguish between + abundances from 
[O II] AA3726+3729 and AA7320+30 because they have very different E ex and n c . 



4.2. O II ORLs temperature. 

The two phase model predicts that T c (ORLs) < T e (He I) < T C (H I) < T c (CELs) (pu] 
2003|) , with the difference between any of these two temperatures being proportional to the 



ADF (see e.g. Figure ||). 

I 1 , . 

Wesson et al. (12003) made use -for the first time- of the temperature sensitive ratio 

J(A4089.29)//(A4649.14) to determine T e of the ionized gas in which O II lines arise in PNe. 
These authors found very low T e in two H-deficient knots in the PN Abell 30. Later, several 
authors have found similar results in other PNe (e.g. [Tsamis et al.| |2004| ; |Liu et aL| |2004 ; 



Wesson et al.|[2005| ). The intensities of the ORLs that originate from states with different 



orbital angular momentum present different dependencies with T e . For example, comparing 
the intensity of one 3<i-4/ transition with respect to a 3s-3p transition, is possible to compute 
the electron temperature. However, this method has its difficulties: first, the dependency 
of the intensity ratio with the temperature is weak, so these lines must be measured very 
accurately; second, the relative intensities of the O II lines can be affected by departures of 
the local thermodynamic equilibrium (LTE) in the fundamental recombination level of the 
0++ ion, 3 P. |Tsamis et ag ([2004]) consider the ratio of the 3d-4/ transition of O II A4089.29 
and the 3p-3s transition O II A4649.14 of multiplet 1 because these lines arise from states of 
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Fig. 6. — ADF (obtained for different ions) vs. excitation energy of the upper level from which 
each line originates (down) and critical density (up). In upper panel, the represented ions 
are, from left to right: 0+ ([O II] AA 3726+29), 0++, 0+ (AA 7319+20+29+30), Ne ++ and 
C ++ ; in lower panel the represented ions are, from left to right: + (AA 7319+20+29+30), 
0++, Ne++, 0+ ([O II] AA 3726+29) and C++. 
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high total angular momentum J, 3 P2 and, therefore they must be affected in a similar way 
by this effect. 

Determination of electron temperature from O II ORLs in H II regions is a non trivial 
problem because in these objects ORLs are usually much weaker than in PNe. In Figure |7| 
we have represented the change of J(A4089.29)//(A4649.14) ratio with T e for T e =10 4 cm 
~ 3 , and we have compared that function with the ratios obtained in NGC 3576 and the 
Orion nebula, the only objects of our sample where O II A4089.29 line has been measured. 
Additionally, we have represented the ratio measured in the giant H II region 30 Doradus 
in the Large Magellanic Cloud (LM C) ( |Peimbert| |2003| ) and in several PNe taken from the 



literature: M 1-42 and M 2-36 ( Liu et al.| |200TD , [Isamis et aT| (|2004l) sample, PNe with 



ADF > 4 of IWesson et aE] (|20Q5|) sample, NGC 5307 ( Ruiz et al.| |2003|) and NGC 5315 



( Pcimbert et al.| [2004|) . This last PN is also included in the PNe sample of [Isamis et al.| 



( j2004|) and we have considered it with the purpose of comparing the effect that high spectral 
resolution observations introduce in the determination of T e (0 II). NGC 5315 is appropiate 
because NLTE effects are small due to their high n e (~10 4 cm -3 ) and also presents an ADF 
which is similar to that measured in H II regions. Within the sample of PNe included in 
panel a) of Figure |7j there are objects with very high ADF (between 5.0 and 22) and in 
panel b) of Figure |7j there are other PNe that show values similar to that obtained in H II 
regions. NGC 5307 and NGC 5315 show a moderate ADF (1.9 and 1.7 respectively), and 
have been studied from high resolution spectra (R~8800), which avoids the overlapping of 
O II multiplet 1 lines (in the other PNe, O II A4649.14 has been measured deblending the 
line using a multiple Gaussian fit) and of O II A4089.29 line with other spectral features^. As 
we can see in panel a) of Figure PNe with high ADFs show, in general, very low T c (0 II), 
which is consistent with the two phase model predictions. On the other hand, PNe with 
moderate ADFs (panel b) of Figure [7]) present temperatures that, in general, are also below 
T c (CELs) and T e (H I), except in the cases of NGC 5882 ( [Isamis et aq P041 ) and NGC 5315 



( Peimbert et al.| |2004|) , where T e (0 II) is consistent with the derived T e from CELs. On 
the other hand, T e (0 II) obtained for NGC 5307 is extremely low, which can be due to 
uncertainties in the measurement of the intensity of O II A4089.29 lineal (see Figure §). In 



4 In some of the PNe of the Tsamis et al. sample, the contribution of Si iv A4088.85 line to the intensity 
of O II A4089.29 line was corrected using the measured intensity of Si iv A4116.10 line and adopting the 
theoretical Si iv A4088.85/Si iv A4116.10 flux ratio (2:1). 

5 The NGC 5307 spectrum was obtained in the same observation run as those of NGC 3576 and the Orion 
nebula. In the analyzed spectra of these observations, O n A4089.29 line is closely blended with a spectral 
feature produced by charge transfer effects in the CCD due to the presence of a very bright line in a different 
order; this feature is clearly separated from the O n A4089.29 line in the Orion nebula, but not in NGC 3576 
nor in 30 Doradus (see Figure @). An alternative form to correct it, would be to adopt the theoretical ratio 
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the case of the Galactic H II regions, T e (0 II) obtained for the Orion nebula is somewhat 
higher to that obtained from CELs; for NGC 3576, T c (0 II) is smaller than that obtained 
from CELs, but O II A4089.29 line is affected by charge transfer effects in the CCD (see 
Figure §). Similarly, J(4089)//(4649)=0.825 ratio found for 30 Doradus (|Peimbert| |2003f ) 
is extremely high, due to a significant contribution of charge transfer to the intensity of 
the line (see Figure U). F° r these two regions we have tried to separate the line from the 
parasite feature, but it was not easy because of the low signal-to-noise of the data in the 
spectral range of interest. As a first step we have identified the line responsible of the charge 
transfer effect in the CCD and then we have measured the flux of all the features produced 
by charge transfer in the CCD (hereinafter CCT) in all the orders; finally, we have fitted 
the flux decrement through the different orders in order to estimate the contribution of the 
CCT to the O II A4089.29 line. Unfortunately, uncertainties in the observed fluxes were 
high, so we can only make a rough estimation to the corrected line flux. For NGC 3576, 
we have estimated that CCT is about a 16% of the observed flux of the O II A4089.29 
line. In Figure |7]c we have represented the new corrected value for NGC 3576, adopting the 
same uncertainty in the line flux that before correction. In the case of 30 Doradus, we have 
estimated that the contribution to the observed flux of O II A4089.29 line is between 30% and 
60%. We have represented an intermediate case (correction of 45%) in Figure [j]d, assuming 
that the uncertainty in the J(4089)//(4649) ratio of 30 Dor is about 30% or even greater. 
As it can be seen in Figure [7|c, applying the estimated corrections, the T e (0 II) obtained for 
NGC 3576 is now consistent with T e derived from CELs. For 30 Doradus, taking into account 
that the n e derived by |Peimbert| (|2003|) was ~300 cm -3 , we have represented the theoretical 
J(A4089.29)/J(A4649.14) ratio for n e = 250 cm -3 (P.J. Storey, private communication) which 
is quantitatively different from that derived for n e =10 4 cm -3 ; unfortunately, in this case, the 
results are not as clear as in the case of NGC 3576: the large uncertainties in the adopted 
flux of the O II A4089.29 line makes the O II 7(4089)// (4649) ratio unreliable to compute 
T e because uncertainties are compatible with variations of T e from 600 K to 25000 K (see 
Figure |7]d). In order to have additional H II region data, we have represented the value 
of the intensity ratio for non-published UVES datasets for two slit positions of the Orion 
nebula: the position labeled as "Orion 1" coincides with one of the two zones studied by 
Esteban et aL| (p98|) , and one slit position located on the Orion bar (24" N and 12" W 
of 8 2 0h A). From the position of the error boxes for these two regions in Figure |7]c, it 
can be seen that the values of T e are similar to that obtained for the Orion 2 slit position. 
Considering T e ([0 III]) calculated for these two regions and the relatively high densities that 
have been computed in both regions -n e ([Cl III]) = 7900±1300 and 5300±800 cm" 3 for 



with other line of the same multiplet but unfortunately we do not know the effect of NLTE to the relative 
intensity of the lines of the multiplet. 
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Orion 1 and Bar, respectively, which minimizes NLTE effects in the intensity of the lines-, 
we can affirm that T e (0 II) measured in different positions from the Orion nebula are similar 
or even somewhat larger than that measured from CELs. On the other hand, it is true that 
in the presence of pure temperature fluctuations, T e (0 II) should be somewhat lower than 
that measured from CELs but, taking into account that measured t 2 's in H II regions are 
moderate, differences as large as those predicted by the two-phase model are not expected. 

Unfortunately, the available data for H II regions are not enough to draw definitive 
general conclusions; however our analysis allows us to enumerate a series of problems that 
we have found and to reach some partial conclusions; a) first, it is necessary to develop 
atomic models that consider departures from LS coupling in NLTE, with the purpose of 
studying how this effect affects each line; b) it is necessary to carry out very deep high 
resolution spectrophotometry in H II regions to extend the catalogue of this type of objects 
with measured T e (0 II); c) the few available data for H II regions behave differently from 
those of most PNe and do not follow the predictions of the two phase model as given in 
Tsamis fc Pequignot] ( |2005| ) for 30 Doradus. 

Finally, a similar diagram has been constructed using C II RLs ratios: /(A9903) / /(A4267) 
and 7(A6462)//(A4267) -that are not affected by NLTE effects- but these ratios are much less 
temperature sensitive than O II A4089.29/A4649.14, and variations between 10 3 to 15xl0 3 K 
are of the order of the observed error, so these ratios are not reliable to determine accurately 
electron temperatures. 

5. Conclusions. 

We have derived the abundance discrepancy factor, ADF, defined as: ++ /H + (ORLs)/0 ++ /H + (CELs) 
in eight Galactic H II regions, finding similar values for all of them. Furthermore, these val- 
ues are similar with the ADF determinations in extragalactic H II regions available in the 
literature. This is a fairly different behavior from the one shown by the ADF in PNe. 

We have found that -within the uncertainties- the ADF is not related to O/H, ++ /H + , 
nor to the ionization degree. As well, it is not related either to the assumed T e (High), so 
systematic effects in the abundance determination from CELs can be discarded. Also, the 
ADF is not related to the T (Low)/T e (High) ratio, suggesting that T e natural gradients in 
H II regions do not produce the abundance discrepancies. 

In the case of H II regions, contrary to the case of PNe, the ADF seems to be metallicity 
independent in the studied range of O abundances. Also, there is no evidence that the ADF 
is related to the effective temperature of ionizing stars -as should be implied by photoelectric 
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Fig. 7. — O II A4089.29/A4649.14 line ratio as a function of the electron temperature. Solid 
line in a), b) and c) panels is the theoretical ratio for n e = 10 4 cm -3 . The data are represented 
as error boxes (see text for the references). In panel a) we have represented a sample of PNe 
with ADF > 5.0; in panel b), we have represented a sample of PNe with moderate ADFs, 
similar to those found in our sample of H II regions; in the panel c), we show the two H II 
regions for which it has been possible to measure O II A4089.29 3d-4f transition: NGC 3576 
and the Orion Nebula and two additional slit positions in the Orion nebula (see text) and, 
finally, in panel d), we show the data for 30 Doradus ( |Peimbert| |2003| ). In this last diagram 
we have also included the theoretical behavior of the line ratio with respect to electron 
temperature for n e = 250 cm -3 (P.J. Storey, private communication) in order to perform a 
more appropiate comparison in the case of 30 Doradus, which shows n c ~ 300 cm -3 . For 
NGC 3576 and 30 Doradus, solid line boxes indicate the values found after correction by 
the presence of charge transfer features in the CCD (see text), the dashed box indicates the 
position of NGC 3576 before the correction. The horizontal arrows indicate the value of the 
line ratio for some objects which corresponding temperature is outside the scale, whereas 
the vertical arrows indicate that the line ratio is outside the represented scale. 
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Fig. 8.— O II A4089.19 line in 5 spectra of H II regions and in PNe NGC 5307. It can be 
noticed that this line is severely blended with another feature produced by charge transfer 
along the CCD (CT) in NGC 3576 and 30 Doradus, whereas in NGC 5307 the signal to 
noise is too low. The Orion nebula is the only H II region in which this line has been well 
measured. The calibrated spectra of 30 Doradus and NGC 5307 have been kindly provided 
by Antonio Peimbert. 
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heating due to the presence of dust near the stars-. We have explored the behavior of the 
ADF with respect to the FWHM of the spectral lines finding no correlation although the 
result is not conclusive due to the paucity of objects with lines of large FWHM. 

We have found that the ADF seems to be slightly related to the excitation energy, E ex 
and to the critical density, n c (in contrast to that observed in PNe), a behavior that is not 
predicted by the two phase model initially proposed for PNe and recently extended to H II 
regions ( [I'samis fc Pequignot| |2005| ), but consistent -at least in the first aspect- with the 



predictions of the standard temperature fluctuations paradigm. On the other hand, electron 
temperatures derived from ORLs are consistent with those derived from CELs, within the 
uncertainties, also contradicting the large differences predicted by the two phase model. All 
these results suggest that the physical mechanism that produces the abundance discrepancy 
in PNe and H II regions might be different. Our results do not permit to discern what is 
the mechanism or natural phenomenon that underlies the abundance discrepancy problem, 
but seem to be more consistent with the predictions of the pure temperature fluctuations 
paradigm than with those of the two phase model. Therefore, it seems necessary to explore 
further other alternatives for the origin of the temperature fluctuations and the abundance 
discrepancy in H II regions. 

Obtaining deep spectra of extragalactic H II regions is fundamental to check the different 
models proposed to explain the abundance discrepancy, permitting to extend the parameters 
space (metallicity, ionization budget, kinematics, morphological complexity) of the objects. 

Finally, we want to stress again that, although the presence and survival of temperature 
fluctuations in gaseous nebulae is a controversial fact, it is clear that the behavior in H II 
regions and PNe should be explored in independent ways, due to their different origin, 
evolution time scales and physical processes that could affect both type of objetcs. 
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Table 9. ADFs, excitation energy and critical density of the CELs used to compute the 

abundances. 



Ion 


log(ADF) a 


E ex (cm x ) 


n c (cm 3 ) 


C++ 


0.35±0.07 


52400 


1.2xl0 9 


Ne++ 


0.28±0.09 


25520 


l.OxlO 7 


0++ 


0.25±0.02 


20000 


6.4xl0 4 


+b 


0.16±0.06 


26800 


4.8xl0 3 


O+c 


0.15±0.06 


13660 


5.8xl0 6 



a Weighted mean for all Galactic and extra- 
galactic H II regions. 

b From [O il] AA 3726+29. 

c From [O n] AA 7319+20+29+30. 



